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SIMPLE  ANALYSIS  OF  LIGHT  ION  BEAM  LOSSES 
IN  DEUTERIUM  PLASMA  CHANNELS 


I.  Introduction 

in  a  previous  paper,1  the  magnetohydrodynamic  response  of  plasma  channels  to  propagating  light 
ion  beams  has  been  investigated.  The  half  a  centimeter  radius  channels  were  assumed  to  be  deuterium 
at  around  10" 5  g/em  with  SO  nsec  pulses  of  ,4-IMA/cm2,  3-5  MeV  proton  beams  injected  into  them. 
Basic  phenomena  like  beam  collisional  energy  losses,  heating  of  the  channel,  expansion  of  the  channel 
due  to  beam  pressure  and  the  resulting  reduction  in  confinement  were  shown  to  take  place  over  the 
duration  of  the  beam  pulse.  While  a  window  on  the  plasma  density  was  given,  no  optimum  for  propa¬ 
gation  was  defined.  In  this  paper,  we  present  a  simple  analysis  to  optimize  power  density  transport  of 
the  ion  beam  through  the  channel  and  show  when  such  an  optimization  process  is  possible,  in  particu¬ 
lar,  we  concentrate  on  the  energy  delivery  at  the  tail  of  the  beam  because  both  voltage  and  current  are 
highest  at  that  time  (as  a  result  of  bunching  requirements.)  We  consider  separately  the  ion  beam  total 
energy  losses  and  the  current  density  reductions  in  the  plasma  channel.  We  look  successively  into 
these  two  areas  and  combining  them  together,  we  look  into  the  optimization  of  the  power  density  after 
transport  in  the  channel. 

This  analysis  applies  primarily  to  deuterium  plasma  channels  where  radiation  losses  are  assumed 
to  be  negligible.  If  the  beam  pulse  length  were  short  enough  or  the  plasma  channel  density  high 
enough,  negligible  motion  would  be  induced  into  the  channel  by  the  beam  over  the  pulse  duration  of 
the  beam  pulse.  However  for  very  short  beam  pulse  length,  not  enough  energy  would  be  available  in 
the  beam.  For  high  channel  density,  beam  energy  losses  would  be  too  large  and  therefore  one  must 
decrease  the  channel  density  and  allow  channel  motion  on  the  time  scale  of  the  beam  pulse  duration 
This  motion  then  changes  the  channel  confinement  characteristics. 
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Two  kinds  of  effects  are  thus  to  be  considered:  energy  losses  and  beam  current  density  reductions. 
The  beam  energy  losses  are  due  to  decelerating  electric  fields  in  addition  to  the  binary  collisional  energy 
losses  and  beam  current  density  reductions  are  related  to  the  plasma  channel  radial  expansion.  These 
two  different  kinds  of  effects  will  be  treated  separately.  Combined  together,  they  lead  to  a  new  concept 
of  overall  transport  efficiency  which  we  define  as 


“  VtVj  ~ 


An 


(1) 


where  J  refers  to  beam  current  density,  «  to  beam  energy  (MeV/proton),  "ouf  to  values  at  the  end  of 
the  channel  and  "in”  at  the  channel  entrance  (beam  injection  values). 


A  brief  examination  of  the  channel  parameters  shows  that  the  plasma  density  in  the  channel  and 
the  plasma  channel  radius  just  prior  to  beam  injection  are  some  of  the  easiest  to  change.  Parameters 
which  characterize  channels  are  chemical  composition,  magnitude  of  the  current  lch  which  flows 
through  them,  current  density  profile  jci,(r),  radius  rcM,  mass  density  distribution  p  (r)  and  also  radial 
velocity  distribution  in  the  channel  v(r)  prior  to  beam  injection.  Recently,  imploding  channels2  have 
also  been  proposed  as  a  medium  for  transporting  ion  beams  and  they  are  presently  under  investigation. 
In  the  present  paper,  we  confine  ourselves  to  the  following  choices:  Deuterium  gas  constitutes  the  back¬ 
ground  channel  material  in  order  to  be  able  to  neglect  radiation  losses.  The  channel  discharge  current 
follows  from  the  beam  characteristics  only  and  is  given  by  /rtC4)  =  .5  10“%£  K0  (cm/sec)  where 9 m  is 
the  maximum  beam  injection  angle  and  Vo  the  ion  beam  initial  velocity.  The  current  density  is 
assumed  to  be  uniform  up  to  the  channel  radius.  In  a  light  ion  beam  transport  scheme,3  4  the  beam  is 
first  focused  then  transported  at  a  radius  comparable  to  the  target  radius  before  hitting  the  target. 
Different  possibilities  for  the  relationship  between  the  channel  and  the  target  radii  are  described  later  in 
this  paper.  The  channel  configuration  has  been  assumed  to  be  a  z-discharge  initially  at  rest  (vr  -  0  at 
t  -  0).  This  has  been  shown  to  be  experimentally  and  theoretically  MHD  stable  on  the  channel  forma¬ 
tion  time  of  a  few  p  sec.5  6  The  mass  density  radial  distribution  in  the  channel  prior  to  beam  injection  is 
assumed  to  be  independent  of  the  radius.  Its  magnitude  can  be  varied  over  several  orders  of  magnitude 
up  to  a  limit  p,  when  energy  losses  become  too  large.  This  value  is  given  by 
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paSL  -  a«  (2) 

where  S  is  the  binary  collisional  stopping  power,  L  is  the  length  over  which  the  beam  is  to  be  tran¬ 
sported,  and  a  is  the  fraction  of  the  initial  beam  energy  which  is  allowed  to  be  lost  during  transport. 
For  fully  ionized  deuterium,  S  =  103/c  MeV/g/cm2  where  «  is  the  proton  energy  in  MeV.  The  tran¬ 
sport  length  L  derived  from  bunching  considerations  is 


where  vo  and  vr  are  respectively  the  ion  velocities  at  the  head  and  the  tail  of  the  beam  and  r  is  the 
beam  pulse  duration,  all  these  quantities  being  defined  at  the  channel  entrance.  For  vr  —  2v0  (energy 
ramp  increasing  to  twice  its  initial  value),  l(cm)  -  3.3  109  «t  (sec).  Combining  Eqs.  (2)  and  (3),  we 
find  a  value  forp0  equal  io 

a  3.10~,3«3/2 

P.  “  - - - •  (4) 

Fore  —  5  MeV,  t  —  50  nsec,  and  a  —  20%,  p„  —  1.3  10“ 5  g/cm3. 

A  minimum  density  must  also  be  achieved  in  the  channel  in  order  to  ensure  charge  and  current 
neutralization  and  also  not  to  excite  electrostatic  instabilities.  If  such  instabilities  are  excited,  the 
plasma  resistivity  becomes  much  larger  than  Spitzer  resistivity  which  causes  an  increase  in  the  net 
current  as  well  as  large  electric  fields.  If  these  instabilities  are  not  excited,  Spitzer  resistivity  can  be 
assumed  and  typical  diffusion  time  of  the  net  current  are  of  the  order  of  a  few  p  sec.  Thus  during  a 
typical  beam  pulse,  a  net  current  of  the  order  of  1%  of  the  beam  current  will  be  added  to  the  channel  as 
was  observed  in  numerical  solutions.1  To  avoid  electrostatic  instabilities,7  np  >  5  x  1017  for  the  above 
parameters.  This  minimum  density  is  about  an  order  of  magnitude  smaller  than  p„  given  in  the  exam¬ 
ple  above,  hence  no  plasma  instabilities  occur  if  more  than  a  few  %  fractional  energy  loss  is  allowed. 

II.  Electric  Field  Generation 

At  (  <  0,  there  is  a  small  electric  field  in  the  channel  corresponding  to  lch.  As  the  beam  enters 
the  channel,  electric  fields  from  various  sources  develop.  From  the  generalized  Ohm’s  law,  they  are 
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due  to  channel  expansion  with  velocity  \p  which  leads  to  vpB/c  and  plasma  resistivity.  In  addition, 
there  is  an  equivalent  electric  field  due  to  beam  energy  coflisional  deposition. 


The  electric  held  terms  are  all  negative  with  the  possible  exception  of  the  first  one  for  an  implod¬ 
ing  channel.  The  energy  losses  can  be  written 


dE 

dx 


V— ‘  +  pS  +  T )jp 


where  the  brackets  mean  average  over  an  ion  betatron  orbit.  The  first  term  goes  from  0  at  the  channel 
center  (B  =»  0,  v  «•  0)  to  a  maximum  value  at  the  channel  boundary  where  vn  and  B  usually  reach  their 
maximum.  For  the  purpose  of  the  present  simple  analysis,  a  mean  value  is  used  for  this  term  which  is 
equal  to  1/2  of  its  maximum  value.  Term  by  term,  the  electric  field  is  thus: 


a) 


From  Ref.  1,  we  get  v 


<jbB 


„  -  .  However,  this  relation  is  correct  to  zeroth  order  only.  If  we  take  into 

Pac 


account  the  fact  that  the  mass  density  decreases  due  to  expansion,  then  to  first  order,  the  expression 
for  the  radial  velocity  becomes 


where 


__  Ob(\B0 
^  cp0(  1+x) 


3  10  2 jbcJihot2 


P  o 


(5) 


(6) 


is  the  expansion  parameter  and  the  0  subscripts  refer  to  channel  and  beam  conditions  at  injection.  The 
reduced  velocity  is  due  to  the  fact  that  the  (j :xi)  term  driving  the  expansion  decreases  with  radius 
faster  than  the  mass  density  (the  current  density  decreases  in  the  same  way  as  the  mass  density  and  the 

if  field  decreases  by  a  factor  (1  -t-  *)  as  will  be  shown  in  Sec.  III).  Now  Bn  “  7^-,  so  that 

5r„ 


<vVg>  _  ~  64  10  " 

f  2c  P0(l+X)>' 


(7) 
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Note  that  this  term  is  inversely  proportional  to  p()  for  small  \  and  proportional  to  p0  for  large  x-  Our 
simple  analysis  does  not  apply  to  that  latter  case. 

b)  pS. 

This  term  is  simply  — — — 7  —  V/cm  for  protons  in  deuterium  with  T,,  >  IeV.  It  is  proportional  to 
(1  +  x)-  « 

p()  for  small  x  and  becomes  proportional  to  p,'  for  large  values  of  x  (which  correspond  to  small  p0). 

c)  nj,,. 

Using  Spitzer  resistivity  justified  by  operating  at  densities  around  10l!t  cm-1  to  avoid  instabilities,  this 
term  is  proportional  to  —'77-  jh  where  the  plasma  return  current  density  jn  has  been  assumed  to  be 

equal  in  magnitude  to  the  beam  current  density  jh.  Due  to  the  small  net  current  (J,h  <  <  jh)  the  most 
critical  quantity  to  determine  in  this  expression  is  Tt..  Assuming  no  radiation  losses,  we  know  from 

Ref.  1  that  Tt.  ~  3.5  10'  so  that  for  x  <<  1. 

«  *  ”  v/cm 

where  /  is  in  sec  and  jw  in  A/cm:.  An  order  of  magnitude  estimate  shows  that  this  resistive  electric 
field  is  equal  to  40  V/cm  for  js,  =  101,  A/crrr,  /  =  50  nsec  and  «  =  3  MeV.  Its  dependence  on  density 
is  very  weak  for  small  x-  Collecting  all  the  terms,  the  total  energy  losses  may  then  be  written  as 

(IE  =  A  Bpo  +  Cy,;  ^ 

dx  p»(i+x(f)l:  U+x('H:  y^:/''-ii  +x<')l 

jr 

From  the  first  two  terms,  we  see  that  —  may  reach  a  minimum  as  a  function  of  p(,  for  small  x-  We 

ax 

find  that  this  minimum  is  reached  for 


1  -X  ■  A 

1  +  3x  B 


Note  that  x  <  1  is  a  necessary  condition  for  a  minimum  to  exist.  For  small  x<  the  density  correspond¬ 


ing  to  this  minimum  is 


p"  "  \B 


u/:  |  4  l0  r 
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This  expression  physically  equates  the  energy  losses  due  to  binary  collisions  to  the  energy  losses  associ¬ 
ated  with  the  plasma  expansion.  Going  back  to  the  general  case,  since  we  have  seen  that  *  is  a  func¬ 
tion  of  p0,  Eq.  (10)  is  an  implicit  equation  for  p0.  Eliminating  this  p0  between  Eqs.  (7)  and  (10),  a 
relationship  between  the  various  beam  parameters  when  a  minimum  in  the  electric  field  occurs  is 

obtained  I  .  _  1 

A0f3-  2.2  10--»  rf2o*x2|_L_^_J  (12) 

This  equation  shows  that  when  a  minimum  E  field  exists,  a  given  fractional  expansion  x  sets  a  relation¬ 
ship  between  the  beam  pulse  duration  and  the  beam  current  density.  A  complete  solution  to  Eq.  (9)  is 
shown  in  Figs.  1  and  2.  It  corresponds  to  jbo  “  1  MA/cm2,  /rt#  “  50  kA,  rcho  —  0.56  cm,  e*  “  3 
MeV  in  Fig.  1  and  fb  —  5  MeV  in  Fig.  2.  We  find  indeed  that  the  minimum  in  the  electric  field  exists 
only  for  small  pulse  durations.  From  Eq.  (12)  it  is  apparent  that  ytgr3  has  a  maximum  as  a  function  of 
X  which  occurs  for  x  ”  6.  For  longer  pulse  durations  or  larger  jb  there  will  be  no  minimum  in  the 
electric  field.  This  happens  for  t  >  25  ns  in  the  example  of  Fig.  1.  At  the  same  time  that  the 
minimum  disappears,  the  expansion  increases  and  while  the  beam  current  density  decreases  the  electric 
field  losses  decrease  also.  For  bunched  cases,  (Jb„t  ™  const),  we  see  from  Eq.  (12)  that  a  minimum 
will  be  reached  more  easily  since  now  / 2  is  smaller.  It  is  important  to  keep  the  resistive  term  in  the 
electric  field  at  shorter  pulse  durations  (unbunched)  since  in  that  case,  the  resistivity  can  be  significant 
because  the  beam  did  not  have  time  to  heat  up  the  channel  to  large  temperatures.  The  expression  for 
the  energy  efficiency  follows  directly  from  Eq.  (9).  It  can  be  written  as 

_  *out  _  *  in  ~  EL  —  j  _  EL 
*  in  *  in  *  in 

Using  typical  £-field  values  of  Fig.  1,  E  —  1.5  kV/cm,  L  -  3m  and  eb  ”  3  MeV,  n,  equals  85%. 

III.  B,  Field  Modifications 

At  t  <  0,  the  B%  -field  configuration  is  assumed  to  be  that  shqwn  in  Fig.  3a.  It  corresponds  to  a 
uniform  channel  current  density  following  the  formation  of  a  low  temperature  channel  where  the  mag¬ 
netic  field  had  time  to  diffuse.  After  beam  injection,  we  assume  that  the  temperature  of  the  channel 
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has  increased  due  to  beam  deposition  and  that  the  expansion  is  self-similar  (the  velocity  increases 
linearly  from  the  center  to  the  edge  of  the  channel).  Because  of  the  large  conductivity  due  to  channel 
heating,  we  assume  that  the  5-field  has  been  convected  out  with  the  plasma  channel  and  that  the  max¬ 
imum  5-field  is  now  reached  at  r(l+x)  as  shown  in  Fig.  3b.  Its  maximum  value  has  decreased  from 
50  to  5^(1  +x)  since  the  same  total  current  fch  is  assumed  to  flow  in  the  channel.  In  fact,  this 
assumption  keeps  the  total  flux  f  Bdr  constant  and  the  beam  stays  confined.  We  assume  that  the  beam 
adjusts  instantly  to  new  channel  conditions  and  that  the  change  in  beam  radius  follows  the  change  in 
channel  radius.  Thus,  because  of  the  expansion,  the  current  density  is  reduced  by  a  factor  (1  +  x)2- 
This  model  assumes  a  free  channel  expansion  and  does  not  take  into  account  the  fact  that  there  can  be 
a  sharp  density  jump  at  the  channel  boundary,  nor  that  the  temperature  has  remained  low  in  the  outer 
region  leading  to  diffusion  of  the  5-field  there. 


The  current  density  reduction  factor  is  thus 


1 


(13) 


A  o  d+x)2' 

In  the  present  form  (no  Vp  effects  on  xh  it  does  not  depend  on  beam  energy.  It  does  not  show  any 
extremum  and  leads  to  low  values  of  Vj  (^50%)  when  channel  expansion  occurs  (x  >  0.4).  This 
current  density  drop  due  to  channel  expansion  may  be  more  restrictive  than  electric  field  losses  as  seen 
in  Fig.  4. 


We  can  now  combine  the  energy  transport  efficiency  v,  and  ny  in  one  quantity  defined  as: 

-  v,Vj- 

This  expression  can  be  written  in  terms  of  first  order  quantities  again  and  the  analysis  is  quite  lengthy. 
It  can  be  shown  that  a  maximum  for  i\„  occurs  as  a  function  of  p0.  Physically,  it  happens  because  v, 
decreases  for  increasing  Po  while  rj y  increases.  An  implicit  equation  for  this  optimum  value  of  the  den¬ 


sity  is  given  by 
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(l-COU+x)2*  ±  V(l-  C02(l+x)V-t-  /f'g'a  +  SyHl  +  ax) 

Po”  /i'd+sx) 

where 

/f’-  AL/e 
B'~  BL/t 

C'  -  CeinL/<Jbot)'n. 

Figures  5  and  6  show  tj„  for  the  same  parameters  as  used  in  Figs.  1  and  2.  We  see  that  the  maximum 
occurs  for  increasing  p0  as  the  beam  pulse  duration  increases  and  that  this  maximum  is  very  broad. 
The  maximum  values  for  i),r  depend  on  the  energy  of  the  beam  which  affects  the  energy  transport 
efficiency.  These  calculations  have  been  made  for  a  3  m  long  channel  —  longer  channels  would  reduce 
the  energy  efficiency  further.  For  these  conditions,  we  see  that  it  is  difficult  to  achieve  an  overall  t),r 
greater  than  50%  for  long  pulses  (t  ^  40  nsec).  In  order  to  increase  the  total  efficiency,  it  is  necessary 
to  consider  transport  schemes  that  have  the  following  options: 

a)  the  channel  expansion  is  limited  by  means  such  as  surrounding  gas  blanket  or  magnetic  field. 

b)  the  channel  implodes  initially. 

c)  the  channel  is  allowed  to  expand  but  is  followed  by  a  focusing  element.8 

IV.  Effects  of  Plasma  Channel  Radius 

The  previous  analysis  has  been  made  supposing  that  the  plasma  channel  radius  was  0.56  cm  ini¬ 
tially  (corresponding  to  an  area  of  1  cm2).  However,  the  conclusions  just  reached  suggest  the  use  of  a 
focusing  element  after  the  transport  plasma  channel  as  one  possibility.  This  solution  has  very  important 
implications  for  the  whole  transport  scheme  concept.  It  indicates  that  a)  constraints  on  the  channel 
radius  can  be  relaxed  and  transport  can  be  performed  at  a  radius  larger  than  the  target  radius  b)  current 
density  reductions  can  be  compensated  for  by  focusing  the  beam  at  the  end  of  the  channel  whereas 
energy  losses  cannot.  In  this  section,  we  look  briefly  at  the  influence  of  channel  radius  on  both  energy 
losses  and  transport  total  efficiency.  The  formulation  is  the  same  as  that  used  previously.  We  just 
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replace  the  beam  and  channel  current  densities  by  the  beam  and  channel  total  current  divided  by  ir  rc2h0. 
The  total  efficiency  can  then  be  written 


n,t- 


1 


1- 


A" _  ,  Bp0  C"rch0 


T  +  " 

T  +  - 

Po^0 

1  +  -T- 

i  +  4- 

1  +  -^ 

>ch0 

o 

■« 

c- 

rch0 

L 


where 


(U) 


«  “  3.10  3  h0lch//P o 
A"-  6.4  x  10 -uVcV 
C"-  1.77  C«,n/V(4l0/2r3/2). 

Note  the  4th  power  dependence  on  r,.*  which  comes  into  various  terms  and  which  has  been  pointed  out 
in  the  expansion  parameter  x  previously.  In  general,  it  is  complicated  to  see  the  effects  of  channel  ra¬ 
dius  on  either  the  energy  losses  or  the  transport  total  efficiency.  Solutions  to  Eq.  (14)  appear  in  Fig.  7 
for  the  energy  loss  efficiency  and  in  Fig.  8  for  the  transport  total  efficiency  as  a  function  of  rthg  and  p0 
for  the  following  parameters.  /„  —  1  MA,  lth  =  50  kA,  e,  —  3  MeV,  /  =  40  nsec,  L  *=  3  m.  The 
effects  of  the  plasma  channel  radius  are  now  more  easily  apparent.  The  energy  loss  efficiency  depends 
weakly  on  the  channel  radius  at  low  channel  density  as  shown  in  Fig.  7.  In  general,  using  small  radii 
causes  the  plasma  to  expand  so  that  the  final  radius  does  not  depend  strongly  on  initial  radius.  The 
curves  drop  as  a  functions  of  p0  because  of  the  collisional  energy  losses. 


As  for  the  overall  transport  efficiency,  note  in  Fig.  8  the  S-shape  of  the  curves.  At  small  radii, 
current  densities  increase  and  channel  expansion  can  be  large  leading  to  small  r At  large  radii, 
current  densities  decrease  as  well  as  channel  expansion  and  y\„  becomes  equal  to  tj,.  Increasing  the 
channel  radius  to  more  than  1.2  cm  at  p0  "  10-6  g/cc  or  to  more  than  0.6  cm  at  p0  -  1.6  10~5  g/cc 
will  not  bring  any  significant  increase  in  the  value  of  rflr.  At  these  radii,  the  decrease  in  y),r  with  p0  is 
due  to  the  decrease  in  v,  with  p0. 
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Notice  that  these  results  are  given  at  r  -  40  nsec  while  previous  results  were  given  as  a  function 
of  time. 

V.  Conclusion 

We  have  shown  that  in  order  to  minimize  electric  field  losses  in  a  plasma  channel  where  a  high 
intensity  light  ion  beam  is  injected  an  optimum  density  exists  for  low-Z  material  and  short  pulses.  We 
have  also  shown  that  no  optimum  density  minimizes  reduction  in  current  density  associated  with  the 
outward  convection  of  the  A-fleld.  Finally,  we  have  seen  that  an  optimum  density  may  exist  for  the 
overall  beam  energy  flux  efficiency. 

A  practical  way  of  choosing  the  initial  channel  density  in  a  channel  initially  at  rest  is  to  ensure  that 
it  remains  close  to  p  „  which  is  the  maximum  allowable  beam  energy  degradation,  thus  providing  the 
highest  inertia  to  beam-induced  motion  and  limiting  changes  in  the  initial  B-fleld. 

In  order  to  increase  the  overall  beam  power  density  transport  efficiency  one  of  the  following 
schemes  needs  to  be  considered:  channel  surrounded  by  gas  blanket,  imploding  channels  and  Anally 
channel  followed  by  focusing  element.  The  advantage  of  this  latter  scheme  is  that  it  can  cancel  com¬ 
pletely  in  principle  the  reduction  in  current  density,  leaving  only  energy  losses  in  the  channel  which 
depend  only  weakly  on  channel  radius.  For  example,  the  case  of  /*  *■  1  MA,  /<•*  —  50  kA,  L  —  3  m, 
r  -  0.56  cm  shows  an  overall  efficiency  less  than  50%  for  beam  pulses  longer  than  40  nsec.  However, 
when  this  same  beam  current  is  allowed  to  propagate  in  a  channel  of  1  cm  radius  for  example  followed 
by  a  focusing  element,  then  the  overall  efficiency  can  get  as  high  as  90%. 

Acknowledgment 

This  work  was  supported  by  the  Department  of  Energy. 


10 


NRL  MEMORANDUM  REPORT  4640 


References 

1.  D.G.  Colombant,  S.A.  Goldstein,  and  D.  Mosher,  Phys.  Rev.  Lett.  45,  1253  (1980). 

2.  S.A.  Goldstein  and  D.A.  Tidman,  Proc.  IEEE  Int.  Conf.  on  Plasma  Science,  Madison,  Wise. 
(1980). 

3.  S.A.  Goldstein,  D.P.  Bacon,  D.  Mosher,  and  G.  Cooperstein,  Proc.  Second  Int.  Conf.  on  High 
Power  Electron  and  Ion  Beam  Res.  and  Tech.,  Ithaca,  N.Y.  (1977). 

4.  P.A.  Miller,  L.  Baker,  J.R.  Freeman,  L.P.  Mix,  J.W.  Poukey,  and  T.P.  Wright,  ibid. 

5.  F.L.  Sandel,  F.C.  Young,  S.J.  Stephanakis,  W.F.  Oliphant,  G.  Cooperstein,  S.A.  Goldstein,  and  D. 
Mosher,  Bull.  Am.  Phys.  Soc.  25,  900  (1980);  J.N.  Olsen,  D.J.  Johnson,  and  R.J.  Leeper,  Appl. 
Phys.  Lett.  36,  808  (1980). 

6.  W.M.  Manheimer,  M.  Lampe,  and  J.P.  Boris.  Phys.  Fluids  16,  1 126  (1973). 

7.  P.F.  Ottinger,  S.A.  Goldstein  and  D.  Mosher,  NRL  Memorandum  Report  4548  (1981). 

8.  P.F.  Ottinger,  S.A.  Goldstein,  and  D.  Mosher,  Proc.  IEEE  Int.  Conf.  on  Plasma  Science,  Madison, 
Wise.  (1980). 


II 


•7 


I)  (i  COIOMBANT 


a)  tsO 


I  i|t  \  -  tal  Magnetic  field  configuration  in  the  channel  helot e  beam  miection 
U  one  \pond  mg  to  uniform  channel  cutfcnt  densin  * 


b)  t:  T 


th)  Magnetic  held  configuration  in  the  channel  afici  beam  inicclion  (corresponding 
lo  same  net  current  and  taking  into  account  channel  expansion* 


14 


7 


Fig  4  —  Current  density  reductions  as  a  function  of  beam  pulse  duration  and  background  channel 
density  for  the  same  parameters  as  in  Fig.  I  < this  curve  does  not  depend  on  ion  beam  energy) 
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Fig  5  —  Total  power  density  transport  efficiency  as  a  function  of  pulse  duration  and 
background  channel  density  for  the  same  parameters  as  in  Fig  I 


unction  of  background  channel  density  and 
ne  parameters  as  in  Fig.  I  at  f  -  40  nsec 
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